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Abstract 

We develop a gauge theory for diffusive and precessional spin dynamics in two-dimensional 
electron gas with disorder. Our approach reveals a direct connections between the absence of the 
equilibrium spin current and strong anisotropy in the spin relaxation: both effects arise if the spin- 
orbit coupling is reduced to a pure gauge SU (2) held. In this case, by a gauge transformation in the 
form of a local SU (2) rotation in the spin subspace the spin-orbit coupling can be removed. The 
resulting spin dynamics is exactly described in terms of two kinetic coefficients: the spin diffusion 
and electron mobility. After the inverse transformation, full diffusive and precessional spin density 
dynamics, including the anisotropic spin relaxation, formation of stable spin structures, and spin 
precession induced by a macroscopic current, is restored. Explicit solutions of the spin evolution 
equations are found for the initially uniform spin density and for stable nonuniform structures. 
Our analysis demonstrates a universal relation between the spin relaxation rate and spin diffusion 
coefficient. 

PACS numbers: 72.25.-b 
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I. INTRODUCTION. 



Description of the spin dynamics of a two-dimensional electron gas is one of the most 
important problems for fundamental and applied modern spintronics.[l-3] Two mutually 
related problems in this field attract a great deal of attention: spin current and spin relax- 
ation. Spin currents describe how the spin density pattern changes with time mainly due 
to the spin transfer between different parts of the electron gas. Since the spin dynamics 
of interest occurs usually in systems out of equilibrium, spin relaxation becomes important 
and contributes strongly into the evolution of spin density pattern. 

The key for understanding these properties is the spin-orbit coupling, making the orbital 
and spin degrees of freedom mutually dependent. Spin-orbit coupling has many crucial 
influences on the properties of the systems where it occurs: the typical examples are nuclei, 
elementary particles, atoms, and electrons in solids. Spin-orbit coupling, in general, makes 
the spin a non-conserved quantity, thus leading to a spin relaxation. It causes mutually 
dependent spin, charge, and mass flows in solids and quantum liquids. ^-6] In addition, 
spin-orbit coupling leads to a spin response to an external electric field, providing an ability 
of spin manipulation by the electric field driving the dynamics in the orbital degrees of 
freedom. j3] 



ib- 
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The general techniques for calculation of spin relaxation and spin current out of equi 
rium are the classical or quantum Boltzmann-like equation for the spin-density matrix 8- 
and nonequilibrium Green functions. 3, Q] In this approach, the description of the electron 
dynamics takes into account possible relaxation processes due to the presence of spin-orbit 
coupling, disorder, phonons, and electron-electron collisions. The experimentally observable 
spin dynamics is due to the spin-orbit (spin-momentum) coupling. At the equilibrium, the 
expectation value of the spin current can be calculated directly. Surprisingly, such a direct 
calculation demonstrated that the spin current can exist even in the equilibrium state of a 
two-dimensional electron gas with spin-orbit coupling. [15] This observation brought a puzzle 
for the understanding of the basic phenomena in spin transport since the equilibrium spin 
current is not related to any spatial spin accumulation that can be seen experimentally. 

On the other hand, the spin dynamics due to the spin-orbit coupling can be understood 
in terms of a theory where the couplin g is treated as a non-Abelian gauge field, and the 
corresponding formalism can be applied \L(L \v\ . On a single-electron scale, for example, for 



2 



electrons in quantum dots, the gauge transformation of the spin-orbit field was employed in 



Refs. 
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19| and used for the analysis of experimental results on spin manipulation by electric 



field in quantum dots in Ref. 20(. Other interesting theoretical examples of applications of 
the non-Abelian gauge field approach for single-electron spin transport and electrons in 
quantum dots, were found and studied (for example, Refs.|21|-|26j|). 

When applied to the two-dimensional electron gas, the approach based on a formal SU (2) 
gauge invariance of the spin-orbit Hamiltonian (i. e. the symmetry with respect to local 
rotations in the spin subspace) proved that the equilibrium spin current is the diamagnetic 



response to the effective non-Abelian spin-orbit magnetic field. [27| If the spin-orbit field is 
a pure gauge and, thus, can be removed by a gauge transformation, the effective SU(2) 
magnetic field is zero, and the equilibrium spin current vanishes. 

Here we present a theory based on the gauge transformation, for spin dynamics in a two- 
dimensional electron gas in the case when the spin-orbit field can be completely removed by 
such a gauge transformation. We show that the absence of the equilibrium spin current is 
directly related to the giant anisotropy in the spin relaxation rate, when the relaxation does 
not occur for certain spin directions. 28|, |29| After gauging away the spin-orbit coupling, 
the entire nonequilibrium dynamics of a transformed spin becomes almost trivial and can 
be described phenomenologically exactly by only two transport coefficients which can be 
determined experimentally, or calculated theretically to any desired level of accuracy. The 
first is the spin diffusion coefficient and the second is the electron mobility required only 
when a constant electric field is applied. With the inverse transformation to the initial 
dynamical variables, we recover the full nontrivial dynamics of the physical spin, including 
the absence of the spin relaxation for certain spin directions, that is a strong anisotropy 
in the spin relaxation, stable spin configurations forming persistent spin helices, and spin 
precession due to a charge current in a constant external electric field. In addition, this 
approach allows making predictions for more general cases of spin-orbit coupling, including 
nonuniform spin-orbit fields. 
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II. SPIN CURRENT AND SPIN RELAXATION: THE CONVENTIONAL AP- 
PROACH. 



We begin with the conventional Hamiltonian of spin-orbit coupling in two-dimensional 
electron gas: 

3 

where aj a is the coordinate-dependent spin-orbit coupling field, kj = —id/dxj is the mo- 
mentum operator, Cartesian subscript indices j — x,y correspond to the electron coordinate 
p = (x, y) , and o a are the Pauli matrices with the upper Cartesian indices corresponding 
to three directions x, y, z in the spin subspace. We use the system of units with h = 1 and 
sum up over repeating indices. Interaction H so arises in a two-dimensional electron gas from 
various sources. Two origins are considered as the most important. The first one, arising due 
to the inversion asymmetry of the crystal unit cell, is described by the Dresselhaus model. 
The other one is the Rashba field, Q where the coupling originates from the macroscopic 
asymmetry of the structure hosting the two-dimensional electron gas. 31] Due to various 
physical origins, including material, structure, doping, and possible mechanical strain, nu- 
merical values of parameters <x, a vary strongly from system to system ranging from 10~ 12 
eV-cm for Si- to 10 -9 eV-cm for GaAs-based structures and will not be discussed here. 

For the coordinate-independent spin-orbit field kjaj a (p) = 0, and the Hamiltonian ([T]) 
can be presented as: 



H so = J2 a i ( h bV °) ( 2 ) 
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where hui is a unit-length vector, and ctj is the corresponding spin-orbit coupling constant 
for given component of momentum; its contribution to the Hamiltonian is, therefore, propor- 
tional to the spin projection onto the hui axis. The coupling leads to a momentum-dependent 
spin splitting of the electron states. As a results, the Fermi line of the electron gas becomes 
spin-dependent and two Fermi lines in the electron gas appear. The H so Hamiltonian makes 
the velocity spin-dependent: 

k- k- 
with m being the electron effective mass. 
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With the spin-dependent velocity in Eq.Q we define the operator of the spin current in 
the form: 



k 

where the SU(2) group generators r a = a a /2, and C£, C*k are the corresponding spinors. 
For example, let |$) be the ground state wave function. The resulting expectation value of 
the total spin current, summed up over all electrons in the gas is: 

< j;) = <*i m ■ (5) 

In the absence of special symmetry relations between the components of the Hamiltonian 
H so , the expectation values of spin current ( Jj 1 ) are not zero, leading to the puzzling equilib- 
rium spin current without any measurable spin transport. Therefore, spin current can be a 
characteristic of the equilibrium states of two-dimensional electron gas. In the conventional 
calculation of (J?) due to the spin-doubling of the Fermi line, contributions to (Jj 1 ) come 
from two subsystems: single- and double occupied states at a given electron momentum. 
These two contributions have opposite signs and almost compensate each other, yielding 
the results in the third order of the coupling constants a|. This third-power dependence 
is expected from perturbation theory: (j") should be an odd function of the spin-orbit 
coupling and vanish in the first order since in the ground state without spin-orbit coupling 
the Fermi-line is not spin-split and all states with given k are doubly occupied. 

Another important feature of the electron gas with spin-orbit coupling is the spin relax- 
ation. Assume that one has initially produced a nonequilibrium state $5 of a uniform spin 
density with the components: 

5 a =($ 5 |^C;r a C k |$ s ). (6) 

k 

Then, the state will relax to the equilibrium through all possible interactions and spin- 
orbit coupling. The first stage of the process, the momentum relaxation, is fast. If the 
spin-orbit coupling is weak compared to the random interactions causing the momentum 
relaxation, as it is assumed for the rest of this paper, the following spin relaxation is relatively 
slow. As a result, at the second stage the spin components decrease with relaxation rates 
described by a symmetric tensor T ab : 

^- = -T ab S b . (7) 
5 



The components of T ab depend on the spin-orbit coupling and all possible interactions of 



electrons with disorder, phonons, and other electrons in the system, 
coupling vanishes, T ab = 0. 
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III. SPIN-ORBIT COUPLING AS A GAUGE FIELD: PURE GAUGE. 



Now we write the Hamiltonian of two-dimensional electron gas in the presence of spin- 
orbit coupling as: 

H= — [ dxd y m + (idi + Ai) 2 m + W (tf + (8) 
2m J 

where W (\I/ + , ^) contains all explicitly spin- independent terms, including electron-electron 
interactions and possibly, the effect of the external potential. The general non-Abelian 
two-component potential is given by 2 x 2 matrices: 



Aj = A*T a = 2ma 3 (p)h a m (p)r a . 



(9) 



The expression (Q is valid for any arbitrary nonuni 
perform at a given p— point a local SU(2)— rotation 



orm spin-orbit coupling. Let us now 
2t| in the spin subspace by 



U = exp[i0 o (p)T°] 5 



(10) 



with the transformation of the field operators: 



(11) 



This transformation renders the spin-independent quantities such as the charge density and 
the charge current density, invariant. In contrast, the spin density operators, 



S = S a r a , 



(12) 



transforming as 



s = usir 1 , 



(13) 



exemplify covariant observable quantities. This difference between the physical quantities 
which transform invariantly and covariantly under a local SU(2) rotation is crucial for the 
understanding of the spin dynamics. 
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For the matrix U = exp [i6 (h • r)] , where h is a unit length vector, the r b — matrices, 
transformed according to Eq. (1131) . acquire the form: 



r b = h b (h-r) + cos #[r b -h b (h- r)] + sin 6e abc h a r c , (14) 

where 

£ abc ig 

the Levi-Civita tensor. This equation shows that the product h r is unaffected 
by the transformation f JTOj) . Therefore, if we present S as the sum of longitudinal and 
transverse components S = S\\+S± with S\\ = h (S ■ h) , the longitudinal component (spin 
projection onto the h— axis) remains constant, while the S± does not; it rotates by the angle 
9 around the h— axis). This simple observation will be important for the further analysis in 
this paper. 

The Hamiltonian preserves its form under a local SU(2) rotation of the fermionic fields 
if the vector-potential is transformed as follows 

Ai = IT 1 (idiU) + IT^U. (15) 

Indeed, after the transformation the Hamiltonian acquires the form: 

H=^-J dxdy^ + (idi + X) ^ * + W (^ + , ^) , (16) 

which is identical to that of Eq. but with \1> and Ai being replaced by the transformed 
quantities, \1/ and Ai, respectively. 

Assume now that Ai in the original Hamiltonian ()8]) corresponds to a pure gauge vector- 
potential, that is both A x and A y can be removed by the above transformation such that 
A x — Ay — 0. In this case there exists a local rotation determined by three coordinate- 
dependent functions 9^(x,y): 

XJ A = expii6 a A (p)T a ], (17) 
such that the initial components Ai can be presented in the form 

Ai, = U A {idi\J A l ) . (18) 

A vector-potential of this form is gauged away by the transformation f lTD]) with U = U_4: 

A = \J- A l {id t + V A (^U^ 1 )) U4 = 0. (19) 

If the spin-orbit field can be removed by a gauge transformation, the subsequent spin dy- 
namics is simplified considerably and in certain regimes, like the drift-diffusion processes 



considered below, the problem becomes elementary. The inverse SU(2)— rotation trans- 
forms the spin components to the actual values and we recover the full dynamics of the 
physical spin. We will follow this procedure in the present paper. 

We mention a textbook example of a similar approach. When the motion of a relativistic 
electron in static perpendicular electric field E and magnetic field H is considered, there 
exists a reference frame, where, after the Lorentz transformation, the smaller of these fields 
vanishes. In this frame the equations of the electron motion are very simple, and in the case 

brmation provides 
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H < E, they are essentially, one- dimensional. The inverse Lorentz trans : 
the full description of the electron motion in the presence of both fields. 

In the pure gauge field after the local SU{2) transformation = exp[i^(p)r a ] the 
Hamiltonian takes the form: 

H = J dxdy^> + M> + W (^ + , , (20) 

with no spin-orbit coupling present. As mentioned above, the spin dynamics with this 
Hamiltonian can be formulated in general terms phenomenologically and then by inverse 
transformation, returned to the form where the coupling and full spin dynamics are restored. 

Vector-potential is a pure gauge, allowing removal six terms in A x ,A y , with the trans- 
formation U^4 based on the three functions 0\{p) given certain relations between the A x 
and A y components. The corresponding conditions are naturally formulated in terms of a 
non-Abelian field strength tensor J-^- : the vector potential is locally a pure gauge if the 
field strength is zero, 

For the spatially uniform case, using Eq.(j9]) this condition is reduced to [Ai, Aj) = 0, that 
is: 

(i) either oiiOij = 0, or 

(ii) [h[j] • r, hy] • r] = if a^j ^ 0. 
The commutation relation 

[h W T, h b1 r] = it ■ (h w x h b] ) , (22) 

demonstrates that the spin projections commute only for the same axis, that is hw = ±hy]. 
Therefore, the solution to Eq. (122]) has the form (we assume below hm = hyi in the case (ii) 
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FIG. 1: (Color online) Illustration of two cases of the pure gauge spin-orbit field, (a) one of the 
coupling constants aj is zero, case (i) (b) both coupling constants are not zero, the directions 
of corresponding magnetic fields coincide, case (ii). The direction of the spin-orbit field remains 
constant for any electron momentum k = (k x , k y ). 

for definiteness) : 

Aj = Imavj (h • r) (23) 
where h = = hyj if a^aj ^ or h = h[/] for nonzero a/, where / = x or / = y, as 

1 /2 

illustrated in Fig.(l). Here a = (a^ + a y) > an d f is a unit vector. The corresponding 
gauge transformation is: 

= exp [2imapji>j (h • r)] exp [2imapiVi (h ■ t)] = exp [2ima (h • r) (p ■ u)] . (24) 

From this condition we immediately conclude that the projection of the total spin at the 
hfji = hyi axis commutes with H so , and, therefore, remains constant with time for arbitrary 
dynamics. Experimentally, this fact corresponds to the vanishing relaxation for this spin 
direction; this conclusion crucial for the design and application of spin-based devices. If 
ajCtj = 0, the problem immediately becomes one-dimensional, trivial from the diamagnetic 
response interpretation of the equilibrium spin current, since one- dimensional systems 
do not demonstrate this kind of response. The same situation occurs in quantum wires, 
where the motion of electrons is strictly one-dimensional, no equilibrium spin current exists, 
and the spin projection along the hrji axis is conserved. In the Appendix, for illustration, 
we perform a conventional calculation of the equilibrium spin current in a two-dimensional 
electron gas with the pure gauge spin-orbit coupling and in a quantum wire, and demonstrate 
that the spin current vanishes in both systems. 

There are two widely studied realizations of the above discussed pure gauge field. The 
CKjOj = case is realized for the Dresselhaus model for the electron gas confined in the 
quantum wells of GaAs grown along the [110] crystal axis. The coupling constant a in this 

9 




FIG. 2: (Color online) Illustration of coordinate-dependent mutual orientation of S and S vectors 
corresponding to Eq. (|26p for a structure grown along the [110] crystal axis. Vectors h and v are 
shown in the Figure. The angle between S and S is determined solely by the x-component of 
p— vector. 



system [35j| is approximately inversely proportional to the square of the quantum well width 
w. In this case the vector-potential and the corresponding transformations are: 

(Ax, A y ) = (2mar z , 0) , XJ A = exp [2imaxr z ] , (25) 

where the z— axis is oriented along the growth direction and the x-axis is that of the unit 
cell. Here we use transformation (|24l) with h =(0,0, 1), v = (1,0), and 9(x,y) = 2max to 
obtain: 

t z = t z , t x = cos 9t x + sin 9r y , r y = cos 9r y - sin 9r x . (26) 

We illustrate the resulting relations between S and S for this simple situation in Fig. (2): 
when S remains constant in space, S turns by the angle 9(x,y) around the z— axis. 

The OLiOLj 7^ case is realized in the compensated Dresselhaus-Rashba model for the GaAs 
structure grown along the [001] crystal axis. Here 

(Ax, A y ) = (2ma (r x - r y ) , 2ma (r x - r y )) , (27) 
U.4 = exp \2ima (x + y) (r x - r y )] . (28) 

Here we obtain with h =(1, -1, 0)/y/2, v= (1,1)/^, and 9 = 2 v / 2ma (x + y): 

r z = cos 9r z - -L sin 9 (r x + r y ) , (29) 
V2 

r x = cos 2 -r x — sin 2 -r v H — sin 9r z , t v = cos 2 -r y — sin 2 -t x H — sin 9r z . 

2 2^/2 2 2 ^/2 

10 



Equations ( 126|) . (|29|) illustrate a general feature of the relations between the original S = 
(S x , S y , S z ) and gauge-transformed S spin densities and vice versa. For the spin-orbit field 
characterized by the direction h, for a uniform coordinate- independent S, the components 
Sy and Sy coincide. The S^-component forms a periodic structure on the spatial scale of the 
order of L so = 1/ma, or h 2 /ma when the units are restored. The mean value (s±(x, y)j = 
for the infinitely large systems considered here, where the boundary conditions do not change 
the spin dynamics. The meaning of the length L so can be understood as follows. Hamiltonian 
([2]) shows that the spin-orbit coupling H so causes for an electron moving with the velocity v, 
spin precession around h with the rate of the order of amv. The corresponding precession 
angle is of the order of amL, where L = vt is the electron displacement. Thefore, L so can 
be viewed as the travel distance at which the electron spin can undergo a full rotation. 
Another circumstance is, however, more important: the spin rotation angle depends only on 
the electron displacement and not on the details of its motion between initial and final points, 
leading to the appearance of stable spin structures, discussed below. Here a numerical value 
of typical L so can be of interest. For GaAs with m = 0.067mo, where mo is the free electron 
mass, and a of the order of 10 -7 meV-cm, L so is of the order of several microns. 

In both these systems, the observed spin relaxation rate is strongly anisotropic with 
one spin component having lifetime orders of magnitude longer than the others. The weak 
relaxation rate for these components is determined by the mechanisms different from the 
homogeneous spin-orbit coupling, most probably, related to the disorder in the spin-orbit 



coupling. [361-138] 



IV. SPIN DYNAMICS: DIFFUSION, PRECESSIONAL BEHAVIOR, AND DRIFT 
CONTRIBUTIONS. 

After the gauge transformation, the spin-orbit interactions is switched off. Therefore, on 
a time scale much longer than the momentum relaxation time, the spin dynamics becomes 
combination of pure spin diffusion and spin drift: 

d t S = DAS + fxEjdjS, (30) 

where D is the spin- diffusion coefficient, \i is the electron mobility, and E is the two- 



dimensional applied electric field 



39] as illustrated in Fig. (3). In Eq. (l3~0~]) we have taken 
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FIG. 3: (Color online) Nonuniform spin density evolving by diffusion and drift dynamics. Small 
square with an arrow illustrates the effect of the external electric field on the nonuniform spin 
dynamics. 

into account that the uniform velocity of electrons is — /xE. These two parameters fully 
describe the drift- diffusive spin dynamics in the absence of spin-orbit coupling. Macroscopic 
motion of electrons (electric currents) can drag nonuniform spin density between different 
parts of the electron gas. This effect leads to the fiEjdjS term in d t S. The initial spin 
density eventually vanishes due to diffusion, however, the total integrated spin polarization 
will remain constant. The diffusive evolution of the transformed spin density DAS occurs 
if the electron free path of the order of I = vr p is much less than the spatial scale of the 
inhomogeneity: i <C L so . This condition can be formulated as fl S oT p <C 1, meaning that the 
spin-orbit coupling is relatively weak. The spatial inhomogeneity of the order of L so and D 
of the order of v 2 r p set the time scale of the diffusion smearing of the S as tp ~ L^ /D on 
the order of Q'^t" 1 , and, therefore, the same spin relaxation time for real spin S. 
The evolution of the physical measurable spin density: 

S = U A SV^\ (31) 

is due to the diffusion and spin precession since the transition of electron from point pi to 
point p2 is accompanied by the rotation of its spin, dependent only on the displacement 
P2 — Pi- Irregular motion in the diffusion process is seen in the spin relaxation, and regular 
drift causes spin precession, with these two processes being mutually related. 

Motion of S is described, therefore, by the following equations for the time evolutions 
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of the spin density, which are obtained by applying the inverse transformation (l3Tj) to the 
drift-diffusion equation fl30|) . 



d t S = D\J A [A (V^SU A )] U^ 1 + (iEjVa [dj (U^SVa)] \J a \ (32) 

The resulting most general equation of motion valid for any pure gauge spin-orbit field takes 
the form 

d t S - DAS - iiEjdjS 

= D{2 [U^VU^, V<S] - 2 (U^VU^ 1 ) S (U^VU^ 1 ) + (U^AU^ 1 ) S + S (AV A ) U^ 1 } 
+ t xE j [(\J A d j V A 1 ),S]. (33) 

The total expression for local evolution of spin density components can be obtained from this 
equation by multiplying both sides by r a and taking the trace using the identity: tr (r a r b ) = 
S ab /2. The result can be presented as the sum: 

d t S a = DAS a + fiEjdjS 11 + BfdjS b - H ab S b - T ab S b . (34) 

The general expressions for non-diagonal tensors of kinetic coefficients entering this equation 
are: 

Bf = -Bf = ADti {r a [Ua^U^ 1 , r b ] } , (35) 

H ab = _ H ba = 2/i £. tr | r a [U^U^, T b ] } , (36) 

T ab = 4D (tr {r a (U^U^ 1 ) r b (U^-U^ 1 ) } - itr {r a (U^AU^ 1 ) r b + r a (AV A ) U^r 6 }) 

(37) 

Now we can study the physical meaning of the obtained non-diagonal tensors and simplify 
the expressions for the time derivatives for uniform spin-orbit coupling with: 

Ua = exp [2ima (h • r) (u ■ p)\ , U^ 1 = = exp \-2ima (hr) (u ■ p)\ , (38) 

With the given form in Eq. fl38|) of and U^ 1 we obtain: 

U^VU^ 1 = -2%ma (h • r) u, A\J A = -4m 2 a 2 V A , AIT 1 = -4m VU^ 1 . (39) 

With formulas ()38l) . (|39l) we obtain for the diffusion-related coefficients: 

Bf = -2m aj De abc h c , (40) 
T ab = Am 2 a 2 D (5 ab - h a h b ) . (41) 
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The corresponding drift-dependent contribution: 

H ab = -mail ■ E) e abc h c , (42) 

describes the spin precession. 

The answer for the spin density S with components (S x , S y , S z ) has the form of three 
terms of different order in a: 

d t S=a t S\ Q +d t 8\ 1 + d&\ 2 . (43) 

These terms have different meaning and can be expressed as: 

d t S\ = DAS + fiEjdjS, (44) 
d t S\ 1 = AmDa {u • V) (h x S) + 2mafi {y • E) (h x S) , (45) 
d t S\ 2 = -4mVD (S - h(h • S)) . (46) 

The d t S\ term describes the standard drift-diffusion spin dynamics for zero spin-orbit cou- 
pling. 

The d t S\ 1 term corresponds to the spin precession due to the spin-orbit coupling. The 

mobility-determined contribution in d t S\ 1 is the precession in the macroscopic spin-orbit 

field arising due to the uniform velocity of electrons. When the electric current is induced, 

the momentum distribution function is shifted such that the momentum has a nonzero value. 

I I 

As a result, the Hamiltonian H so forms a macroscopic spin-orbit Zeeman field [40:] and, as 
a result, a regular spin precession dtS =2ma[i {y • E) (h x S) . If [y • E) = 0, contributions 
of the momentum changes along the x and y-axes in the macroscopic spin-orbit "magnetic" 
field compensate each other, and no regular precession occurs. Thus, Eq. (l4"5l) reproduces 
the diffusive and non-diffusive spin precession. 



The d t S\ 2 term is the Dyakonov-Perel' mechanism of spin relaxation, 41( which can be 



seen from the fact that D is determined by (v 2 ) t p , where v is the electron velocity (see, also 



in Ref. 13j). Taking into account that electron momentum is mv, one can see that <9jS| 2 
corresponds to the Dyakonov-Perel relaxation with the relaxation rate on the order of a 2 k 2 r p . 
The obtained relation between the spin relaxation rate and diffusion coefficient is universal. 
For two different systems with the same sample-dependent ma parameter, the ratio of T ab /D 
remains constant. Since the parameters T ab and D can be measured independently, this 
universality can be verified experimentally. For example, in the measurements performed at 
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the same sample at different temperatures, the ratio T ab /D is expected to remain constant 
for degenerated and non-degenerated electron gas. 

Eqs. (j4"5]) and (j4"5]) show that Sy = h • S does not change with time, as expected, and the 
entire dynamics is solely due to the S^— component. As a simple illustration we consider 
the evolution of an initially homogeneous spin density. By solving equations (I43p - (j46p with 
the initial condition S(p, t = 0) = So we find the spin dynamics 

S(t) = h(h • S ) + {cos(f^t)[S - h(h ■ S )] + sin(fi B t)(h x S )} e- rt , (47) 

where Qe — 2amfi(u ■ E) is the precession frequency in a drift-induced spin-orbit Zeeman 
field, and T = 4a 2 m 2 D is the diffusion related relaxation rate. From Eq. ( 147|) we see that the 
spin precesses with the frequency Qe about the h-axis and its transverse component decays 
at the rate T in such a way that the projection of the spin at h remains stationary. By 
comparing the characteristic time scales of the drift-induced precession and the diffusion- 
induced relaxation, we can estimate the external electric field at which the role of the 
drift-dependent terms becomes important; in particular, the precession becomes visible at 
the scale of the relaxation time. From the condition of still visible precession Qe ~ T we 
find the corresponding electric field E ~ amD/fi. In this field, the precession rate QIe is of 
the order of Q 2 q t p , making the contributions of regular and random motion in the precession 
angle of the same order. If the spin diffusion is dominated by the impurity scattering, 
then D and /i are proportional to the momentum relaxation time r p , and this electric field 
is disorder-independent. However, it can change with the temperature since in the non- 
degenerated gas D approaches the electron diffusion coefficient {sqJ and, therefore, by the 
Einstein relation D = fiT. 

Another interesting effect of a spin-orbit coupling, which follows straightforwardly from 
our formulation - the existence of stable spatially inhomogeneous spin configuration. It is 
easy to verify that a general stationary (<9 t S = 0) solution to the equations (T4^1) -(l46l) is of 
the form 

S(p) = h(h • S ) + cos(2ma(i/ • p))[S - h(h • S )] - sm(2ma(u • p))(h x S ), (48) 

where So is an arbitrary constant vector. This spatially inhomogeneous stationary solution to 
the drift-diffusion equation arises due to the symmetry of the system. As it was demonstrated 
for the particular case of the model with balanced Rashba and Dresselhaus couplings, the 
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symmetry can protect electron spins from relaxation [29J and allows for the persistent spin 



helix structures 



42 



43] of the form of Eq. fj48|) . The fact that the shape of this configuration 



does not depend on the diffusion coefficient D shows that the persistent spin structure is 



insensitive to the spin-independent disorder, in agreement with Ref . [42| The analysis of the 
spin helix stability in the presence of disordered spin-orbit coupling can be found in Ref. 44( It 
is interesting to note that the helix structure is also insensitive to the presence of the electric 
field and, therefore, to the mobility and presence of a transport charge current, at least in the 
linear Ohm's law regime. This seemingly counterintuitive result follows from the fact that 
the drift of the helix governed by the second term in (14"4"|) is exactly compensated by the spin 
precession in the current-induced effective Zeeman field, the second term in (1451) . A similar 
cancellation occurs in the diffusion channel. A diffusive spreadout of the helix, the first 
term in (PHI) , and the relaxation of the transverse component of the spin, ( pl6|) . are balanced 
by the "gradient-precession" contribution, the first term in ( I4"5j) . The persistent spin helix 
configuration ( l48l) has an extremely simple interpretation in terms of the transformed spin 
density S. The general stationary solution of the standard drift-diffusion equation ( )30l) is 
simply a constant S = So- The relation between the physical and transformed spin density 
components yield the conservation So ■ h = S • h. The perpendicular So,± is transformed 
according to Eq. (!3"T!) as S = U^iSU^ 1 with = exp [lima (h • r) (p ■ v)\ from Eq. (|24|) 
according to Eq. ( JT4|) . The sum of the transformed terms is precisely the persistent spin 
helix of Eq. (HE}. 

It is also instructive to look at the precession and relaxation of a spatially homogeneous 
spin S(t), Eq. (HTl) . from the point of view of the dynamics of the transformed spin density 
S. The initial condition S(p,t = 0) = So for the physical spin is mapped to the initial 
configuration for S in a form of a spin helix that is similar to Eq. (PJHJ). The subsequent 
evolution of S is governed by the standard drift-diffusion equation (!30|) . Therefore the 
dynamical behavior is obvious -the initial helix for the transformed spin moves with the 
drift velocity v drift = ma\i{y ■ E), and washes out diffusively. When transformed to the 
physical spin, the drift of the helix translates to the precession, while its diffusive decay is 
mapped to the relaxation of the physical spin. This interpretation clearly explains why the 
relaxation time of the transverse components of the spin is universally determined by the 
diffusion coefficient. The relaxation of the physical spin components is gauge-equivalent to 
a purely diffusive process of washing out the initial helix configuration of the transformed 
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spin. 



V. CONCLUSIONS. 

We developed a gauge theory of macroscopic spin dynamics in a two-dimensional electron 
gas when the spin-orbit coupling can be described as a pure gauge, and, therefore, removed 
by a local SU (2) rotation in the spin subspace. We have shown that for a pure spin gauge, 
equilibrium spin current vanishes and a selected axis of conserved spin projection appears 
simultaneously, demonstrating gauge-related symmetry relation of these effects. After re- 
moving the spin-orbit coupling, we considered macroscopic phenomenological equations of 
spin dynamics, including spin diffusion and spin drift in an external electric field. By the 
inverse SU(2) rotation we obtained the full system of partial differential equations for the 
time- and spatial measured spin density evolution. This system reproduces the physics of 
spin precession, stable spin configurations such as the persistent spin helix, and the result- 
ing strongly anisotropic spin relaxation. Since we described the system without spin-orbit 
coupling phenomenologically, our approach is valid at any temperature and electron concen- 
tration. It predicts that the ratio of the spin relaxation rate to the spin diffusion coefficient 
remains temperature- and electron concentration-independent if the coupling constants do 
not depend on these two system parameters. 

We presented explicit equations for the spatially uniform spin-orbit coupling and their 
solutions describing stable nonuniform structures, the precession and the relaxation of uni- 
form spin polarization. These equations can be explicitly generalized for nonuniform two- 
dimensional electron gas in macroscopic systems. We mention two of them. The first one 
is the GaAs quantum well grown along the [110] direction with a modulated width w(x), 
where the spin-orbit a(x) originated from the Dresselhaus coupling, varies as l/w 2 (x). The 
corresponding spin-orbit field A x = 2m«(i)r 2 , A y = 0, with dA x /dy = remains a pure 
gauge. The other example is the balanced Rashba-Dresselhaus model with the coordinate- 
dependent Rashba and Dresselhaus parameters remaining exactly equal or exactly opposite 
everywhere. As in the [110] structure, variation in the Dresselhaus term is due to the con- 
trolled variation in the structure width, while the control of the Rashba coupling is achieved 
by a coordinate-dependent bias across the well. A different kind of inhomogeneity occurs 
in mesoscopic systems where the effect of the boundary conditions for the coupled spin- 
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charge dynamics becomes important. 451447] Generalization of the gauge theory approach 



for the dynamics at the sample boundaries can be an interesting extension of our analysis 
for the infinite systems. Spin dynamics in these systems is of interest for the fundamental 
understanding of spin transport and for applications in spintronics devices. 
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VII. APPENDIX. 



Here we show by a conventional calculation of the spin current that it vanishes at the 
equilibrium in the considered above pure gauge spin-orbit coupling in a two-dimensional 
electron gas and, similarly, in one-dimensional quantum wires. We begin with the pure- 
gauge two-dimensional Hamiltonian, 

h = H + ax (/iva) kx+ hi +ay (/iva) kyi (49) 

where h is unit length vector and a x , a y are the corresponding spin-orbit coupling constants. 
The spectrum of electrons described by Eq. (j4"9|) is the sum of k x and /^-dependent terms for 
the two spin branches " + " and 

k 2 + k 2 

e± (k x , k y ) = v ± {a x k x + a y k y ) . (50) 

Eq. ( l4*9l) demonstrates that the system with pure gauge spin-orbit coupling remains in a 
certain sense, one dimensional: spin-orbit coupling does not couple different components 
of momentum in the spectrum. For illustration we consider only the x— component of 
momentum and velocity: 

= - ± a Xl v x = t[H,x} = ^ + a x (h a a a ) , (51) 
ok r m m 
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; '+" branch 



"-" branch 




FIG. 4: (Color online) Scheme of the spin-orbit split states in a quantum wire. Ep is the Fermi 
energy. 



yielding the spin current component: 

k 

Taking into account that {cx a ,cr 6 } = 2§ ab , we obtain the anticommutator: 
2 k J 2 ^ m J m 



(52) 



(53) 



The total spin current is the sum of contributions of two subsystems (J^) = +(«/£)_ = 

2 (Jx) + ■ Taking into account that for given branch (& b ) ± = ±h b /h, the (Jx) + spin current 
component becomes: 

de (k) 



2 / " ^ / r7/,- r " 



(54) 



where the integration is perfromed over the area in momentum space occupied by electrons 
from the branch. The value of this integral is zero since this area is restricted by the line of 
the constant Fermi energy Ep. 

For one- dimensional case the situation is the same. We take the Hamiltonian: 



H= + a (h a a a ) k. 

2m 

The eigenstates of this Hamiltonian form two branches: 

k 2 



2m 



± ak, 



(55) 



(56) 



corresponding to two parabolas with the minima at — ko and ko = am, respectively, as shown 
in Fig. (4). To calculate the spin current directly, we perform integration over momenta and 
summation over spin branches. The ground state expectation value is: 



-k +k f 



ko—kp 



m 





rko+kp i 


h b ^j dk + 






J ko—kp > 
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k 



m 



a") +h b )dk 



(57) 



where kp is the Fermi momentum determined by the total concentration of electrons n as 
kp = irn/2. We obtain 



J b ) = 2k F h b + 



2m 



k—dk — 
h 



2 h b 

2k F h b -krh 

m h 



—kQ—kp 
pK . 



ko+kp i 6 

k—dk 

■k +kp h 



(58) 



The minimum position k = am, yields ( J b ) = 0, as expected. 

The absence of the spin current in a wire can be related to the vanishing persistent 
Aharonov-Bohm spin current 



48 



49( in a ring with spin-orbit coupling as the ring radius 



goes to infinity. Indeed, if a ring is pierced by a small-radius solenoid with a finite magnetic 
field flux confined inside it, the field at the ring is a pure gauge, and an Aharonov-Bohm 
current in the ring arises. The spin-orbit coupling in a ring can be understood in terms 
of a spin-dependent gauge, and a corresponding spin current can be induced. This spin 
current vanishes in the R — > oo limit. In general, however, these two objects have a different 
topology: a ring can demonstrate a diamagnetic response, while a wire cannot. 
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